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Abstract 

Integrability conditions for systems of bosons or fermions with seniority conserving hamil- 
tonians are derived. The conditions are shown to be invariant under a large class of trans- 
formations of the interaction matrix elements. Previously published integrable models are 
shown to satisfy these conditions and the existence of a new class of integrable models is 
demonstrated. The number of free parameters in the interaction in these models equals the 
number of single particle levels plus 3 . Equations for the energy eigenstates and eigenvalues 
are derived and the eigenvalues of the complete set of two-body integrals of the motion are 
given for the new class. Some two-body correlations in these eigenstates are derived from 
the integrals of the motion. 

I. INTRODUCTION 

Shortly after the development of the BCS theory of sup er conduct ivity[|^], the Hamil- 
tonian and the associated computational techniques were adapted to create the pairing 
model of nuclear structure. This model has been successful in qualitatively accounting for 
nuclear properties such as the energy gap in the single-particle spectrum of even-even nu- 
clei, odd-even mass differences, and the moments of inertia of deformed nuclei in spite of 
the fact that the methods used by BCS do not apply to finite systems since they violate 
number conservation. A great deal of effort was invested into attempts to restore number 
conservation to the theory. At about the same time, an exact analytical expression!^, |^ for 
the eigenstates and eigenvalues of the pairing Hamiltonian with constant interaction ma- 
trix elements was published and largely ignored by the nuclear physics community. Shortly 
thereafter the corresponding many-boson problem was also solved^. More recently, the 
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attention of the condensed matter community has been directed towards understanding 
pairing correlations in finite systems by a series of experiments on metalhc grains with a size 
of a few nanometers ^ and the exact solution has been extensively utilized in studies of 
these systems^. A separate and important recent development has been the demonstration 
that the pairing Hamiltonian, with constant interaction matrix elements, is an integrable 
system. []^ This system is then both solvable and integrable. Another important result is 
the demonstration of a class of solvable and integrable models involving both pairing forces 
and an interaction term quadratic in the particle number operators[§ ^. This paper draws 
much of its inspiration from these works and extends their results. 

The structure of this paper is as follows. We start in Sec. 2 with a definition of the 
class of Hamiltonians that we will consider. They are the most general, seniority-conserving 
Hamiltonians with two-body forces, which describe systems of fermions or bosons. These 
Hamiltonians have a diagonal single-particle kinetic energy plus a pairing term and a term 
quadratic in the particle number operators with arbitrary matrix elements in all terms. 
We then consider a set of operators that have the same structure as the Hamiltonian and 
derive a set of conditions on the matrix elements in the members of the set that must be 
satisfied if these operators are required to commute with the Hamiltonian. We obtain an 
explicit solution for the matrix elements in the two-body part of these constants of the 
motion in terms of the matrix elements in the Hamiltonian and a set of conditions which 
the Hamiltonian matrix elements must satisfy. Solutions of these conditions yield integrable 
Hamiltonians. We show that, due to a symmetry of the conditions, each particular solution 
gives rise to a large family of solutions. The number of free parameters in each family 
is greater than the number of single-particle energy levels. We also show that previously 
published ^ solvable and integrable models satisfy our conditions. By generating a new 
particular solution of the conditions, we generate a whole new family of integrable models. 
In Sec. 3, we show that members of this new family of integrable models are also solvable by 
giving analytical forms for all their eigenstates and eigenvalues. In Sec. 4, we calculate the 
occupations of the single-particle levels and some two-body correlation functions in these 
eigenstates. 
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II. HAMILTONIAN AND CONSTANTS OF THE MOTION 



We consider the class of seniority conserving Hamiltonians for Fermi or Bose systems 
that can be written in the following form H — T + U + V , where 



r = J]2Q (2.1) 

i 

U = ^ Uij (Qj + 2 s rii) (Qj + 2 s nj) 

In these expressions, i labels a single-particle energy-level of energy and g is an inter- 
action strength. The interaction matrices u and v are real and symmetric. The statistics 
factor s = — 1 (+1) for Fermi (Bose) statistics. The operators n and h are defined by 
Fermions: 

The degeneracy of level i is 2f2j with the factor 2 coming from time-reversal degeneracy 
and a — ±1, • • •, ±0j with 



^ 2 £ + (2-2) 

a=l 

bi ^ ^ (li,—a^i,a 
a=l 

Bosons: 

There are 2 possibihties for bosons. In the first case, the single-particle states are not 
time-reversal cigenstates and the above definition applies. In the second case, the single- 
particle states are time-reversal eigenstates as in the k = state or the bound states of a 
potential. We then define the operators 



— - t 



bi y= ^ ^ (^i,a^i,a 

^ a=l 



(2.3) 
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The single-particle operators a satisfy (Fermi)Bose (anti)commutation relations 
which can be used to derive the commutation rules for the n and b operators 



(2.4) 



The effects of degeneracy are expressed explicitly in the presence of the Qj in the commu- 
tation relations and the effects of statistics in the presence of the s in U and in the factors 
{flj + 2s rij) in the Hamiltonian and the commutation relations. The seniority operators are 
defined by 



and they commute with the Hamiltonian. They count the number of unpaired particles 
in each single particle level. The form of the interaction in (2.1) is chosen so as to mimic 
that of a Heisenberg XXZ spin chain model after a quasispin substitution for the fermion 
pair operators. 

The eigenstates of such Hamiltonians factor into a product of a state of the un- 
paired particles times a state of the paired particles. The energy of the state is the sum 
of the energies of the unpaired particles, which interact only through U, plus that of the 
paired particles. The influence of the unpaired particles on the paired particles comes about 
through through a modification of the degeneracies Qj - blocking in the case of fermions and 
antiblocking for bosons. In what follows, we only consider the seniority-zero paired particles 
and assume that the degeneracies are effective degeneracies which depend upon which state 
of the system is being considered. 

Following Ref. 7, we look for two-body constants of the motion K that have a form 
similar to that of the Hamiltonian, K — X + Y + Z , with 




a=l 
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X = J]2x,n,- (2.5) 

3 

where y and 2; are real and symmetric matrices to be determined in terms oi x , e , u , 
and V so that K commutes with the Hamiltonian. Since x is a set of free parameters, equal 
in number to the number of degrees of freedom of the system, the solution of the resulting 
equations will provide a complete set of constants of the motion specified by a choice of x 
and it will be shown that these constants, with different choices for x , commute with each 
other. 

The calculation of the commutator [H , K] follows from (2.1), (2.4), and (2.5) and 
some tedious algebra. The result is 



[H, K] 

= 9^{-{Ck- 6) Vk,i + (cfc - ei) Zk,i} (bl hi - h\ 



(2.6) 



+ ^ [vk,i {yj,k - yj,i) - {uj,k - Uj^i) Zk,i + {vj^kZj,i - Vj,iZj,k)] bl {Qj + 2s rij) hi - h\ {Vtj + 2s rij) bu 



j,k,l 

where 



^k = Xk - - s g yk,k 



e-k^^k- -j^s g Uk,k 



(2.7) 



Setting the commutator equal to zero yields the equations 



Zk,i = 7 ^Vk,l , kytl 

(cfc - ei) 



(2.8) 



and 
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Vk,i {yj,k - yj,i) - {uj,k - Uj,i) Zk,i + {vj,kZj,i - Vj^iZj^k) = , kj^l (2.9) 

We assume that v ^ and first consider the case in which u = . Then (2.9), with 
j ^ A; , is 



Vk,i {yk,k - Zk,k - yk,i) + Vk,kZk,i = , k^l (2.10) 
Without loss of generality, we can require 

yk,k^Zk,k^O (2.11) 

Then 



yk,i = ^-^^^v,,i ,k^l (2.12) 
[^k — ei) 

which requires that Vk,k be a constant which can be taken to be one. Equation (2.9), with 
j , k and I distinct, and, using (2.8) and (2.12), becomes 



which requires vj^k — 1 and is the simple pairing Hamiltonian with constant interaction 
matrix elements. This is the result of Ref. 7. 

We now assume that both u and v . We can solve (2.9), which is antisym- 
metric in k and I , for y by setting j — k which gives 



Vk,i {yk,k - Zk,k) - Vk,iyk,i + {uk,i + Vk,k - Uk,k) Zk,i ^0 , ky^l (2.13) 
Without loss of generality, we can require 

yk,k = Zk,k = Vk,k - Uk,k = (2.14) 
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Then, using (2.8), we have 



yk,i = ^^^^u,,i ,k^l (2.15) 

Equation (2.9), with j , k and / distinct, and using (2.8) and (2.15) provides an 
integrabihty condition on u and v . This condition is hnear in the x's and, equating their 
coefficients to zero, leads to the single equation 



{x-y)u {x, z) V {y, z)-{x- z)u {x, y)v {y, z) + {y - z)v {x,y)v {x, z) ^ , x y z ^ x 

(2.16) 

where we have introduced the notation 



Uj^k^u{ej,ek) (2.17) 
Vj,k = V {ej,ek) 

and write {x,y,z) for {ej,ek,ei) . Note that this equation is independent of particle 
statistics and single-particle degeneracies. 

Eq. (2.16) has a curious symmetry that allows the generation of a whole class of 
solutions from any particular solution. Assume that we have a particular solution of (2.16), 
u {x, y) and v (x, y) , then there will exist a class of solutions generated from this particular 
solution of the form 



u' {x, y) = J^^ZJ^f (^) ' / (^)] (2.18) 

-'(^>^) = 7(|^7(^W/(^),/(y)] 

where / is some well behaved function of its argument. This result follows from the fact 
that the form of Eq. (2.16) is invariant under the transformation 



u^u' , v^v, x^ f (2.19) 
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Note that iteration of the transformation (2.18) is equivalent to a redefinition of the 
function / so that each class of solutions is closed under this operation. If we require the 
functions / to be monotonic with single-valued inverses, then the transformations have a 
group structure. In the following we will look for particular solutions of (2.16) which will 
provide the basis for classes of solutions generated by the transformations (2.18). The 
members of each class are specified by the set of numbers fi = f (cj) so that the number 
of free parameters in a class is equal to the number of single-particle energy levels plus the 
number of free parameters in the particular solution. 

We first show that recently published[|, ^ results are solutions of (2.16). The first 
class of solutions is based upon the simplest solution of (2.16): u = v = 1 . With this solu- 
tion, f/ is a function of the total number of particles, and the Hamiltonian is equivalent to 
the simple pairing model Hamiltonian with constant interaction matrix elements. The mem- 
bers of this class, after the transformation (2.18), are the descendents of this Hamiltonian 
with 



X — y 

u{x,y) =v{x,y) = —— — - (2.20) 

f{x)-f {y) 

This class can also be obtained directly from (2.16) by taking u = v and dividing (2.16) 

by V (x, y) v (x, z) v {y, z) to get the equation 



{x-y) {x- z) ^ {y- z) _^ 
v{x,y) v{x,z) v{y,z) 

which is linear in 1/v . The solutions of this equation are {x — y) /v (x, y) = f (x) — / {y) 
which is (2.20). A second class of solutions is obtained by writing 



u{x,y) = v{x,y)w{x,y) 

and obtaining from (2.16) 



(x — y),.(x — z),.(y — z) 

i^^^ X, z) - ^ {w X, y) + ^-^ ( = 2.22 

v{x,y) v{x,z) v{y,z) 

Solutions of this equation areH, 0] 
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w (x, y) = Cosh [q {x - y)\ 

q{x-y) Cosh [q {x - y)] 
u {x, y) = - 



Sinh [q (x — y)] 

Note that the first class of solutions is just the q ^ hmit of the second. 
A third and new class of solutions is obtained by writing 



v{x,y)^^{x)^{y) (2.24) 

where 7^ {x) is a polynomial in x and u {x, y) is a symmetric polynomial in x and y . 
Substituting this form into (2.16) gives the equation 



{x — y)u {x, z) — {x — z)u (x, y) + {y — z) 7^ (x) = (2.25) 
which is linear in u and 7^ . The solution of this equation is 



(x, y) = 1 + 7i (x + y) + 72 X y (2.26) 

7 (x) = (1 + 2 71 a; + 72 x^) ^ 

with 7i and 72 free parameters. Note that before the transformation (2.19), the first 
two terms in u can be ignored since they result in terms in the Hamiltonian that are either 
constants, proportional to the total number of particles, or a modification of the single- 
particle spectrum. However, they are important after the transformation. Note also that 
this solution reverts back to the first one for the case 72 = 7i • The case in which there is 
no constant term in u can be handled by making the replacements g — > 5170 , 7i ^ 7i/7o , 
i = 1, 2 , and then taking the limit 70 — > . The net effect of which is the replacement of 1 
in (2.26) by . 

Further study of Eq. (2.16) has not produced any other simple solutions that are not 
the result of symmetry transformations apphed to the solutions (2.20), (2.23), and (2.26). In 
the next sections, we consider the eigenstates of the integrable model (2.26) after a symmetry 
transformation . 
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III. EIGENSTATES 



It is a semiempirical theorem that the seniority-zero eigenstates of integrable Hamil- 
tonians of the form (2.1), have the form of a product of correlated pairs. This is true for the 
states of the Hamiltonians with interaction matrix elements given by (2.20) or (2.23) [|], |^ 
and we now show, by construction, that it is also true for the solution (2.26). 

Assuming that the eigenstates can be taken as a product of correlated pairs, we 
write, for a typical eigenstate of 2N particles, 

N 

\^P) = X{Bl\0) (3.1) 

r=l 

where 

Bl = Y.^r{e,)h] (3.2) 

3 

Then, commuting the Hamiltonian past the product of 's , we have 

N \ N / N \ 

H~2Y,E,-E^'A\i,) = Y.\ n Bl\{[H,Bl]-2E,Bl]\Q) (3.3) 

+ 5 E ( n St)[[ff.i3t].B.]|0) 

p,g=i(P7^'?) \''=i,(»'7^p,g) / 
where we have written the energy of the vacuum state as 



^^°^ = -t5Z«m^^^^- (3-4) 

and where the parameters 2Ep are pair-energies to be determined. The energy of the 
state will be given in terms of these parameters. The commutators can be calculated from 
(2.4) and (3.2) with the resuhs 

{[iJ,i?t] -2E,i?t}|0)=2^ 

3 

(3.5) 



u 



'3,k 



— En 



+ 
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[ [H , Bl] , Bl] = -gs {uj,k [(pp {ej) (pq (e^) + (pp (cfe) (pg {ej)] - Vj^k [(pp {ej) (pq {ej) + (pp (cfe) (pq (cfc)]} b]bl 

(3.6) 

For the new set of integrable models, we take u and v to be given by (2.26). Then, 
after the symmetry transformation (2.18), 



Vj,k = (ej, efc) = 7 (e^-) 7 (e^) ^ — ^7 (/j) 7 (/fc) 



(3.7) 



(/, - fk) 



0p (Cfc) = y 0p (/^) ^ ^^''^ 



with a;p a parameter to be determined. For (3.5), we then have 



{[H,Bl]-2EpBl}\0)^2Y^ 



(/, - A) 



(3.8) 



+fE ^^e, - e,) ^fi^^p (/,) 6] 



|0) 



^ (^p, fk) _ ^ 
{fk - i^p) ^ 



<Pp{e,)b]\0) 



and for (3.6), after algebraic simplification. 



[ [H , Bl] , Bl] = -gs '!}'^'_^^'\ ^ (e, - e,) [0, (/,) 0, (/,) - 0, (/,) 0, (/,)] b]bl (3.9) 



Substituting these results into (3.3) yields 
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N 

N / N 

n 

p=l \r=l,{rjtp) 



(3.10) 



^3 



( fi. — a;„) (a;r, — u>a) 



^ [Jk- ^p) 



|0) 



Setting this equal to zero yields the equations 



+ 9" E 



2 (fk — ^v) (<^o — <^») 



0, p=l,---,Ar 



(3.11) 



^p — —-^ "feCfc 



2^ (fk-^p) 



The energy of the state is then given by 



(3.12) 



N 



E^E^'^ + 2j2Ep 



p=i 

N 



(3.13) 



= £<»)-9EE""=»#^ 

p=i k ^•l'' P> 

For fk — , this expression simplifies to 

Af TV' N 

E = + 2 ^ - t/s ^ (u;p,u;^) - ^1 ^1 ^^fe^i (c^p, e^) (3.14) 

p=l P,9=1(97^P) P=l 

The equations (3.11) can be transformed into a form that allows a two-dimensional 
electrostatic analogy for fermions with s = — 1 Some algebraic manipulation of (3.11) yields 
the equations 



2 + 5[0 + 2s(iV-l)](7i +72u;p) 1 



2 V (A - ^p) 



AT' 

^ E 

g=l(q^p) 



[UJn 



Up) 



0, p=l, 



(3.15) 
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For fermions, these equations can be interpreted as the equations for the unstable equi- 
hbrium positions of unit hne charges at the locations Up in the complex plane interacting 
with charges —^7^/2 at the locations fk and two charges located at the zeros of 7^ {uo) . 
In this analogy, the statistical repulsion between fermions is mapped onto the electrostatic 
repulsion of like charges in two dimensions. Such equations have been shown to lead to the 
BCS equations plus corrections in a carefully defined N 00 limit. |]TD[ Eq. (3.14) then 



shows that the energy of the states has terms proportional to the dipole and quadrupole 
moments of the free charge distribution. For bosons, we have a statistical attraction be- 
tween particles and the electrostatic analogy does not work. Such equations have been 
shown to yield a generalized Bose condensation into the lowest two single-particle levels for 
a repulsive interaction in the N —>■ 00 limit. The Bogoliubov approximation, modified to 
accommodate this generalized condensation, then gives good results in this limit. 

The same techniques can be used to calculate the eigenvalues of the integrals K . 
We consider the solution (2.26) after a symmetry transformation. We then have 



K 

where 



(3.16) 



If we make the choice Xk ^ 5ik and K ^ Ki , then Ki and its eigenvalues, ki , are 



(3.17) 



Ki = 2ni + g ^ 



ifi — fk 



{fi, fk) (bl bk + bi 6i) - '^u ifi, fk) (a + 2 s m) {Qk + 2snk) 

(3.18) 



ki = 9 



N 



^ (/i - ^p) 



= -'-n. J2 Qk 



k{k^i) 



U ifi, fk) 
ifi — fk) 



This result will be used in the next section to calculate the occupations of the single- 
particle levels and some two-body correlation functions in these eigenstates. 
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IV. OCCUPATIONS AND CORRELATIONS 



We can use the eigenvalues of the Hamiltonian and the K ,s and the Hellman- 
Feynman|]l3 theorem to calculate certain two-body correlations in the eigenstates presented 



in the previous section. The Hamiltonian and constants of the motion are functions of the 
parameters in the potential, g , 'ji , and 72 as well as the single-particle spectrum ej , and 
its transform fj = f (cj) . Since the function / is unspecified, we can treat ej and fj as 
independent parameters. In this picture the Up are functions of g and the fj , through 
(3.11), and do not depend upon the Cj . Furthermore, comparing (3.13) with (3.18) yields 
the result 



h - K^'^ = ^.{E- (4.1) 

In order to calculate the mean occupations of the single-particle levels, we use (3.18) 
to write 



n^ = \iK^-9^-^] (4.2) 



2 V dg 

The Hellman-Feynman theorem then yields 



{Mn^\^l^) = \{k,-g^^ (4.3) 



The Rp satisfy the linear set of equations 



N' 



Rp-2s Y 7^^^^. = 1 (4-4) 

77-V {Uj.p — UJq) 



obtained by differentiation of (3.11) with respect to g . 

The Hellman-Feynman theorem can also be used to calculate the expectation values 
of the following operators 
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Ki,j = [u {h fj) {bl bj + b] k) - ^7 ifi) 7 ifj) {^i + 2s n,) + 2 s n,)} 



7(/.) 



or 



K. 



i',3 



2^ 



The expectation value of dKi/dfj is given by 



dK, 



dkj 



and 



(4.5) 



(4.6) 



(4.7) 



ifi fj I 

2^ 



7 (/j) 7 (/i) ^fej 

7(/0 9/.- 7(/.)9/, 



= -n,%7 (/,) 7 (/,) 
where we have written 



p=l 



+ 



R 



■p,t 



(4. 



du). 



(4.9) 



and where the Rpj satisfy the hnear set of equations (4.4), with the right-hand-side 
replaced by 7^ {cUp) / {fj — ujpf , obtained by differentiating (3.11), 



V. CONCLUSIONS 

We have derived a general set of conditions for the integrability of the most general 
seniority conserving Hamiltonian with two-body forces for bosons or fermions. They have 
been used to obtain a new class of integrable models. We have shown that the eigenvalues 
and eigenstates of the Hamiltonian and the associated integrals of the motion can be obtained 



15 



by the solution of a fairly simple set of coupled algebraic equations. While this solvability 
depends upon the integrability conditions in a rather complicated way, we have not been 
able to construct a general result that integrability implies solvability. 

The results of this paper can be generalized to other systems such as a system of nu- 



cleons with an isospin-independent J = pairing in jj-coupling|]TT| or a spin- and isospin- 
independent L = pairing in Is-coupling. []T^ Another system to which these methods can be 
applied is the system of fermions with an equal mix of pairing in the ^5*0 and ^Pq channels. lO 
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